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SUMMARY 

An approximate linearized solution is presented for the wave 
drag and lift of an airfoil generated by rotating a thin supersonic 
profile about an axis parallel, or nearly parallel, to its chord. 

The aerodynamic coefficients are obtained from a surface distribu- 
tion of sources, of strengths proportional to the local airfoil 
slopes, about a cylinder whose radii’s and chord equal those of the 
original ring airfoil. This source distribution satisfies the 
bo unda ry conditions when the part of the wing within the forward 
Mach oone from a point on the wing surface departs only slightly 
from a plane. The solution is therefore accurate for ring airfoils 
having chords that are small in comparison to the radius of rotation. 

The lift coefficient of thin supersonic ring airfoils, based 
on the airfoil-surface area, is one-half the Ackeret value for a 
two-dimensional wing of infinite span. The drag coefficient is 
equal to the sum of the Ackeret value for the given profile (with 
the ring airfoil at zero angle of attack) and the induced drag 
coefficient. These coefficients are probably within 5 percent of 
the correct linearized values for ring airfoils whose chord-radius 
parameter (chord divided by the product of the radius and the 
cotangent of the Mach angle) is within the range from 0 to 0.20. 


INTRODUCTION 

The linearized solution for the aerodynamic coefficients of a 
thin supersonic airfoil whose chord elements form a flat surface is 
currently in an advanced stage. Conical superposition has been 
extensively used to determine performance when the plan form is 
bounded by straight-line segments. Surface source distributions 
have been used to develop analytical and numerical methods of solu- 
tion (references 1 and 2) for an arbitrarily shaped plan-form 
boundary. However, airfoils whose chord surface is not flat have 
been comparatively neglected. 
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Wave drag and lift are evaluated herein for the supersonic 
ring airfoil generated by rotating a thin sharp- edged profile 
seotlon about an axis parallel, or nearly parallel, to its chord. 
The chord-radius ratio of the ring airfoil Is assumed to be small. 
A surface source distribution, as originally proposed in refer- 
ence 3 for essentially flat airfoils. Is used as the basis for the 
analysis. The limiting chord-radius ratios for which this source 
distribution gives a valid linearized solution are estimated by 
comparing the results for the external surface of the ring airfoil 
to those obtained by the numerical method of reference 4. The 
ring airfoil in subsonic flight is considered in reference 5. 

The Investigation, conducted at the NACA Cleveland laboratory, 
was completed during December, 1947. 


SYMBOLS 


The following symbols are used in this report: 


°D 

Of 


'P 


drag coefficient based on wing surface, ( , \ 

^ (Jptj2 Z*ro) 

skin-friction drag coefficient based on wing- surface area, 

( friction drag\ 

\ pU 2 2jtrc / 


lift coefficient based on wing-surface area. 


pressure coefficient, ^u^mta^pressurej 


Jp " 8 


wong chord 


D drag 

L lift 


M free-stream Mach number 

£ local source strength per unit area 

r mean radius of airfoil 
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a area of integration 


t maximum profile thickness 


t/o 

IT 

V 



thickness ratio 

free-stream velocity, taken in positive x direction 
perturbation velocity normal to source surface 

Cartesian coordinate system 



cylindrical coordinate system 



a. angle of attack, radians 

Oq angle between chord of profile section and axis of ring 
airfoil, radians 

3 cotangent of Mach angle, (Vm 2 -i) 

X local slope of wing surface with respeot to direction of 
free-stream velocity 

p free-stream density 

a local slope of wing surface with respect to chord 

cp perturbation-velocity potential 

Subscripts: 

e external surface of ring airfoil 

1 internal surface of ring airfoil 
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METHOD OF AHALYSIS 

The linearized equation for the perturbation- velocity potential 
of an irrotational, compressible fluid is given by 



A general solution of equation (1) can be obtained from a surface 
distribution of sources. The potential at a point in the flow field 
is the sum of the contributions of the elemental sources in the 
forward Mach cone from that point. If the source surface is in the 
x,y plane, the potential at a point in the x,y plane is given by 



where q is the function defining the source- strength distribution. 
The coordinate system and the limits of integration are illustrated 
in figure 1(a) . 

The induced normal perturbation velocity at a point on, or an 
infinitesimal distance from, a plans containing sources is simply 
a function of the local source strength and equals qjt. (See 
reference 3.) Thus, If a planar distribution of sources is to 
represent the flow about the top or bottom surface of a thin super- 
sonic airfoil whose chord surface is flat, the boundary conditions 
are satisfied by setting 



(3) 


where X * v/U is the local slope of the wing surface In idle free- 
streem direction. If the sources are assumed to be in the x,y plane, 
the equation for the velocity potential becomes 
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The following sign convention for X will he used: When the surface 
normal projecting Into the stream has a component in the positive 
x direction, the slope is considered negative; when the surface 
normal has a component In the negative x direction, the slope is con- 
sidered positive. For a flat plate at angle of attack a, the slopes 
for the upper and lower surfaces are constant and have the values 
X« -a and X= a, respectively. 

The potential flow about a thin-ring airfoil can he found from 
a source distribution about a oircular cylinder whose axis is par- 
allel to the free-stream direction and whose radius and ohord equal 
those of the original airfoil; however, the function defining the 
source- strength distribution remains to be determined. In the 
limiting case of a thin-ring airfoil having a small chord-radius 
ratio, the wing section in the forward Mach cone from a point on the 
wing departs only slightly from a plane. From the previous dis- 
cussion, a source distribution defined by the local wing slope 
q. = Mj/tc adequately satisfies the boundary conditions. This rep- 
resentation is directly analogous to the solution of an essentially 
flat wing in that allowance is made for variations of wing slope, 
but the departure of the wing suxfaoe from a plane is considered 
negligible. The velocity potential for points on a ring airfoil of 
small chord, expressed in cylindrical coordinates, is then 



where rdoo, cor, and Sr have replaced dr), rj, and y of equa- 
tion (4). The coordinate system is illustrated in figure 1(b). 

The normal perturbation velocity at a point on a three- 
dimensional source suxfaoe is not a simple function of the local 
source strength. Equation (3), and consequently equation (5), 
therefore do not satisfy the boundary conditions when the chord- 
radius ratio of a ring airfoil is sufficiently large that the 
deviation from a plane of the chord surface in the forward Mach 
cone is no longer negligible. 


WAVE ERAS AND LIFT OF gn^rtTPRICAL-BING- AIRFOIL 

v. 

A cylindrical-ring airfoil is a ring airfoil of zero profile 
thiokness, camber, and flare (oq = 0). The velocity potential 



6 


MCA TW Wo. 1678 


of such an airfoil may be found if the function defining X at all 
points is known, for small angles of attack and a small chord, 
this function is shown in appendix A to be 

- a sin to (6) 

» a sin (0 (6a) 


The expression for the velocity potential of the external surface 
of a cylindrical-ring airfoil results from the substitution of 
equation (6) into equation (5): 


^e 



(-a sin to) &|d&? 

- > 2 <* - "> 2 


rx 


Ua 

JC 


it 


re+£i L 

P? 


0 - 


pr 


sin tod Ct) 




p 2 (0 - CO) 5 


( 7 ) 


It is desirable to solve equation (7) direotly for &p e /dx. in 
appendix B, tills solution shows that 


3*6 13a, sin 0 T 1 /_x\ 2 _1_ /_x\ 4 

"S? " p A “ 4[0rJ + 64 


e 


( 8 ) 


The pressure coefficient for a point on the external surface, based 
on the linearized Bernoulli equation, is then 


C 


P>© 


1^2 
u dx 


-2ct sin 0 

P 


1 




(9) 
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The pressure coefficient for the internal surface is the negative 
of that for the external surface. 

The lift and drag coefficients, referred to the airfoil- surface 
area, may he obtained from the respective summations of the effec- 
tive lift drag forces acting on each element of area: 



The results for the lift and drag coefficients show that 



The series in equations (9), (12), and (13) converge very 
rapidly for values of and c/(pr) that are sufficiently small 

to justify the use of a source distribution defined by the looal 
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wing slopes. No appreciable error will be introduced by assuming 
these series to equal 1. The pressure coefficient then becomes 
simply the Ackeret -value corresponding to the local slope: 

C P,e “ | sln ®) < 14 > 

This result indicates that the spanwise variation of slope in the 
forward Maoh cone, due to chord-surface curvature, may be neglected 
when calculating the pressure coeff ioient at points on an airfoil 
by means of a source distribution defined by q = XU/n. Therefore 
equation (14) could have been obtained directly from equation (7) 
by using the mean value sin 6 instead of sin ( 0 . 

Similarly, the lift and drag coefficients become 

°L"f (15) 

and 

C D = y (16) 

Both coefficients are one-half the Ackeret values for a flat-plate 
airfoil of infinite span. 


WAVE DRAG AND LIFT OF BING AIRFOIL OF GENERAL PROFILE 

In the general case, a ring-airfoil profile seotion has camber 
and thickness and its chord makes a finite angle with the axis of 
the airfoil. The local airfoil slopes are then given by the 
expressions 


X e = °e _ (°o + a sia (17) 

^i “ a l + (oq + a, sin ©) (17a) 

For a small chord, the pressure coefficient has been shown to equal 
approximately the Ackeret value corresponding to the local slope. 

c p,e = | [°e “ (ocq + a sin 6)] (18) 

Cp,i = | [<*. + (ccq + a sin 0)] 


(18a) 
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The lift coefficient obtained by substituting these pressure 
coefficients in equation (10) equals that of the cylindrical-ring 
airfoil. Thus the lift coefficient of a ring airfoil is independent 
of its profile slope. 

The drag coefficient is 


°D - f 2 + ^ +£j o * ai 2 ) ax (19) 

The first term of equation (19) defines the induced drag due to 
lift and is one-half the Aokeret value for induoed drag coefficient. 
The last two terms represent the drag of the ring airfoil at zero 
angle of attack and equal the drag coefficient of an Ackeret airfoil 
(having the given profile section) at angle of attack <xq. 


LIFT-ERAG- RATIOS OF RING AIRFOIL 

A calculation can be made to illustrate the order of magnitude 
of lift-drag ratios obtainable from a ring-type airfoil. For a 
symmetrical diamond profile, no flare, and a skin-friction drag 
coefficient C f , the lift-drag ratio is given by 



( 20 ) 


For small values of t/c and Cj, the L/D performance approaches 
that of an Ackeret airfoil. With increasing t/c, the l/d 
approaches one-half the value for the corresponding Aokeret airfoil. 
Equation (20) can be maximized to give the angle of attack at which 
(L/Dj -mqy occurs for a given t/c and PCf. The result shows 


a “ V 2 (o) + 2 Cf 

and the corresponding value of (L/D) max is 


( 21 ) 
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Equation (22) is plotted in figure 2. The values of alxfoil lift 
coefficient at which (l‘A , ) ma x oocur s are cross-plotted In this 
figure. As shown by equations (15), (21), and (22), this lift 
coefficient is the reciprocal of (l/d)^. 

Little information is available on the value of friction drag 
coefficients at supersonic speeds. The mean value Of <= 0.006 is 
sometimes used for thin supersonic airfoils. If is assumed 

independent cf Mach number, figure 2 indicates that ( L A > ) maz 
decreases vlth increased flight speed, hut this effect becomes 
less pronounced with larger values of t/c. 

Because the flare was considered zero and a diamond profile 
has the lowest possible wave drag for a given thickness ratio, the 
values for (L/Djjaax presented in figure 2 are the largest values 
theoretically obtainable from a ring-type airfoil. 

VALUES OF c/(PF) FOE WHICH AERODYNAMIC 
COEFFICIENTS ARE ACCURATE 

The values of chord- radius ratio beyond which equations (15) 
and (19) are no longer valid can be estimated by comparing the 
external lift and drag coefficients as indicated by these equations 
to the values found from the numerical method of reference 4, which 
is an extension cf the solution for slender, pointed-nose bodies of 
revolution presented in references 6 and 7. A line distribution of 
sources and doublets is placed along the axis of an open-nosed body 
cf revolution, and the strength distribution is so adjusted that 
the flow and cross-flow components of the free-stream velocity 
follow the external contour of the body. The resulting values for 
external wave drag and lift may be considered the correct linearized 
solution. The internal flow, however, is not properly described. 

By use of the line-doublet distribution of reference 4, the 
lift coefficient for the external surface of a cylindrical-ring 
airfoil, expressed in the form 0 /a, was found to be solely a 
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function of the airfoil parameter o/(fSr). The relation is shown 
in figure 3. The decrease in lift with increase in chord indicates 
that the downstream sections of the cylindrical-ring airfoil have 
progressively lower external lift. Equation (15), "based on the 
surface source distribution, yields a constant value for PCi^ e /a, 
but the discrepancy between the two methods is less than 5 percent 
for values of o/(pr) up to 0.20. 

Computations for the external lift and drag of finite-thickness 
ring airfoils indicate that the effect of thickness ratio on the 
discrepancy between the two methods is small for the magnitudes of 
thickness ratio permitted by the linearized theory. These discrep- 
ancies depend primarily on the value of o/(0r). This result is to 
be expected, as ohord-plane curvature, and not thickness ratio, is 
the source of the difference. If the discrepancy for the internal 
surface is assumed to be of the same order as that for the external 
surface, the generalization may be made that the coefficients 
expressed in equations (15) and (19) are correct to within 5 per- 
cent for values of c/(gr) up to 0.20. The degree of error may be 
estimated by the plot of figure 3. 


SUMMARY OF RESULTS 

An approximate linearized solution for the aerodynamic coef- 
ficients of ring airfoils having small chord- radius ratios indicates 
that the lift coefficient, based on the wing-surface area, is inde- 
pendent of the profile shape and equals one-half the Ackeret value 
for a two-dimensional wing of infinite span. The drag coefficient 
equals the Ackeret value for the given profile (with the ring air- 
foil at zero angle of attack) plus the induced drag coefficient. 
These coefficients are probably within 5 percent of the correct 
linearized solution for values of chord-radius ratio from 0 to 0.20. 


Flight Propulsion Researoh Laboratory, 

national Advisory Committee for Aeronautics, 
Cleveland, Ohio, April 30, 1948. 
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APPENDIX A 

FUNCTION DEFUSING LOCAL SLOPES ON CYLINDRICAL RING AIRFOIL 

A cylindrical ring airfoil of small chord and at a small angle 
of attack is assumed. 



for which 

At vertical displacement between corresponding points on leading 
and trailing edges 

Ar radial distance between corresponding points on leading and 
trailing edges 

Then 

Av = c sin a 
Ar = At sin go 

a c sin a sin go 

But 

X = — = sin a sin go 
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Inaamucb. as a is small 


X = a sin. a) 


and with the proper sign convention 


o X s* -a sin co 

JJ ® 

O 

and 


X^ = a sin co 
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APPENDIX B 


SOLUTION FOR ^3? FROM TELOCITY- POTENTIAL EQUATION 
dx 


From equation (7) 


/~x 


cp „ Ua 
e - 


re- 


4| 




sin codto 


Je- 


x-g 

Pr 




(Bl) 


p 2 (e - to) 2 


Because sin to Is independent of £, the order of Integration 
of equation (Bl) can "be reversed and the equation integrated with 
respect to |. The resulting expression for cp e can then he dif- 
ferentiated with respect to x to yield a line integral equation 
for 5cp Q /5x. Reference 1 shows that the result is equivalent to a 
line integration along the leading edge 1=0 of the airfoil. 

Therefore, 


5fe _ Tfa 
the « 


re+x/C pr) 


sin todto 


iM - p2 (e ■ 


Da 

*0 


re+x/(pr) 


sin todco 


Je-x/(pr) 

Integrating hy parts gives 


l /-“ 2 + - e 2 ~ JlL 


(B2) 


_ D a 

•** 

[it sin 6 cos + 

'■'6+x/(pr) 

008 “ sin £7^] da> " 

H p * 



L J 

e-*/(P?) 
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Let 


e - a? _ 

x/(ffr) 


Then 

r-e+x/Or) 


6-x/(pr) 


cos co sin 


■i K- M 1 

|_x/(pr)J 


dco = - 


x 

Pr 


J 


-1 


cos fe 


(•■») 


i'" 1 

x 'cos 6 cos + sin e sin — r j sin" 


PS? 


r (• 


Pr 


n ^-1 

Pr7 


And by sine and -cosine series expansions 


- JL cos 0 f" 1 - _i/“ f + -L /«\ 4 - . . . 

Pr -'l j_ 2 ‘\pr/ 4 * \Pr/ 

x 4 a f" 1 1/ zx ^ 1 f zx\ 3 1 . . 

Pr \ pry 3J\pr/ j 


But 


r- 


sin - ^- zdz = 0 


for n equal to a positive even integer or zero, and 


sin -1 zdz 


zdz 

(B4) 


1 zdz 


(B5) 
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f Z n sin" 1 zdz er ( -\ + X ^ x — 4 . . , i\ 

J 1 n+l\ n + l n-l A n-3 2/ 

for n equal to a positive odd integer. Therefore 

pe+x/(pr) r 


cos co sin 


0-x/(pr) 


— aco = — £. sin e ~ (-1 + — ^ 

x/(pr)J &r PrVS/^ 2/ 


- ~ 2 f i + 5 x l\ _1_ /x_\ 5 « / _ 5 x 5 x l \ 

31 \^pr J 4x 2/ + 5 i 6 \‘ 1 + 6x4x2/‘ 


jr sin 0 cos 1 + 

Pr 


JC.N _1 „ 1 /_x\ _1_ 5X1 

&) 21 X 2 _ 4! X 4 X 2 


./ X \ 6 1 w 5 X 3 X 1 

\pr J Gl 6 X 4X2 " * * * 


Substituting equation (B6) in equation (B3) gives 

d^e TJ a . x . x 

bx * P pr pr 


1 . /x 

\ 4 1 .3X1 

2 ^ pr 

) 4! X 4 X 2 

xf 1 

„ 5 X 3 X 1 

JSr) 6i 

6X4X2 

/ x \4 

/ — Y x 

■** 

1 

■ ■■ , • • • 

Ve*; 

64 
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Chord-radius ratio, c/(p?) 


Figure 3. - External lift coefficients of oylindrical-ring 
airfoils as indicated by line-doublet distribution of 
reference 4 and by surface source distribution. 
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